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We investigate symmetry properties of the solar magnetic field in tfie framework of the 
I— { ' hnear Parker's migratory dynamo model. We show that the problem can be mapped onto 

• the well-known quantum-mechanical double-well problem. Using the WKB approxima- 

tion, we construct analytical solutions for the dipole and quadrupole configurations of the 
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generated magnetic field. Our asymptotic analysis within the equatorial region indicates 
the existence of an additional weak dynamo wave which violates the Hale's polarity law. 
We estimate the spatial decay decrement of this wave. We also calculate explicitly the 



> 

^ splitting of the eigenvalues corresponding to the dipole and quadrupole modes. The real 

' part of the dipole eigenvalue is shown to exceed the quadrupole one. A sufficiently long 

time after generation the dipole mode should dominate the quadrupole one. The relevant 
observational evidences of the properties obtained are discussed. 



o 

p^ . 1 Introduction 



2 . The problem of solar magnetic activity is an old and intriguing subject. A 

\ great amount of observational data have been accumulated since the first 

' half of 17th century. A lot of efforts have been made to explain theoretically 

these remarkable data. However, the problem still bears many unanswered 
^ , questions. Among them, the problem of symmetry of the solar magnetic 

^ ' field exists. This problem is visualized as a symmetry violation in some 

tracers of the solar activity at low latitudes in certain periods of times such 
as the Maunder Minumum (Harvey, 1992; Ribes and Nesme-Ribes, 1993). 

It is clear now that solar magnetic field generation is connected with a 
turbulent motion of the differentially rotating electrically conducting solar 
plasma. The full set of the corresponding magnetohydro dynamic equations 
is too complicated to handle analytically. However, a significant simplifi- 
cation arises if we limit ourselves to considering the large-scale magnetic 
structure only. The equations describing the large-scale magnetic field in a 
thin convective shell were obtained by Parker on phenomenological grounds 
back in 1955. These equations have been intensively investigated for almost 
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fifty years. A more solid basis of describing mean magnetic field was pro- 
vided by Steenbeck, Krause and Radler (1966) who formally derived the 
equations of mean- field magnetohydrodynamics. The Parker equations fol- 
low from the general equations of mean-field magnetohydrodynamics being, 
however, only the leading approximation with respect to the short length 
of a dynamo wave. For a more consistent treatment of the dynamo prob- 
lem, the next-to-leading approximation is required and the Parker equations 
must be slightly modified (see e.g. the recent paper of Galitski and Sokoloff, 
1999). 

There is no doubt that the large-scale magnetic field observed on the 
today's Sun corresponds to the steady regime. Recent investigations of the 
Parker equations in the non-linear case have revealed the following three 
features of the steadily oscillatory solutions (see, e.g., Meunier et ai, 1997; 
Bassom et al., 1999; Kuzanyan and Sokoloff, 1996). Firstly, the structure 
of the steady state solution is not very sensitive to the explicit form of the 
nonlinearity introduced into the mean-field equations. Secondly, the spatial 
profile of the dynamo wave in the steady regime may retain some qualitative 
features of the one for the kinematic problem. Thirdly, the frequency of 
the magnetic field oscillations coincides in the main approximation with 
the imaginary part of the eigenvalue of the linear problem. Thus, it is 
quite reasonable to start up with a relatively simple linear case in which 
the equations may allow for an analytical solution. Moreover, the kinematic 
problem is a first step to approach the nonlinear case. 

In the kinematic approximation, the mean-field equations reduce to an 
eigenvalue problem for a linear differential operator. Formally, this opera- 
tor is quite similar to the standard Hamiltonian in non-relativistic quantum 
mechanics. The diffusion terms play the role of kinetic energy, while the 
alpha-effect (helicity coefficient) and the differential rotation play the role 
of a two-component external potential. It has been found fruitful to ap- 
ply some well-established quantum-mechanical approaches to the problem 
at hand (e.g., Sokoloff et al, 1983, Ruzmaikin et al, 1990). Using the 
WKB approximation, Kuzanyan and Sokoloff (1995) and later Galitski and 
Sokoloff (1998, 1999) have obtained asymptotic analytical solutions of the 
dynamo problem in the framework of the Parker model. The solution de- 
scribing the wave of magnetic activity, the so-called dynamo wave, was built 
up of the three parts (dynamo waves): a strong dynamo wave propagating 
towards the equator in the main spatial region, a relatively weak dynamo 
wave propagating in the subpolar region polewards and an extremely weak 
dynamo wave reflected from the pole and decaying exponentially propagat- 
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ing equatorwards. The former two waves are known from the observations 
of the Sun (Makarov and Sivaraman, 1983), while the latter wave predicted 
theoretically has not yet been observed probably due to its weakness. 

A major assumption adopted in the works cited above was that the 
sources of magnetic field generation in the southern and northern hemi- 
spheres were well separated and the generation took place in different hemi- 
spheres absolutely independently. Certainly, this assumption is not appro- 
priate when considering the symmetry problem in which the interplay be- 
tween the dynamo waves in diflcerent hemispheres is the key factor. More- 
over, the WKB approximation breaks down in the very vicinity of the equa- 
tor as shown below. 

In the present paper, wc investigate the symmetry problem by solving 
the dynamo equations in the Parker's model. We investigate the equatorial 
region and show that the wave of magnetic activity generated in a given 
hemisphere does not vanish at the equator but rolls over it propagating 
ahead in the other hemisphere and decaying with propagation. Due to this 
fact, there is a weak interaction between the dynamo waves in different hemi- 
spheres which results in a small splitting of the eigenvalues corresponding to 
different symmetry configurations of the global magnetic field. We estimate 
the spatial decay rate of this straying dynamo wave. We also calculate the 
splitting explicitly and find that the growth rate for the dipole configura- 
tion exceeds the one for the quadrupole configuration in the framework of 
the Parker's model. Notice, that the asymptotic method used enables us to 
estimate just the asymptotic order of magnitude of this eigenvalue splitting 
and we can hardly calculate the exact numerical coefficient by this method. 



2 Basic Equations 

In the present section we describe the model and establish the necessary 
notations. Since, the current paper is based heavily on our previous work, 
we highlight the main result only briefly omitting all intermediate calcula- 
tions. We refer an interested reader to the original paper of Kuzanyan and 
Sokoloff (1995) [see also Galitski and Sokoloff (1999)] for a more exhaustive 
presentation. 

The starting point is the mean-field equations derived by Krause and 
Radler (1980): 

— = Vx(aB) + Vx(vxB)+/3AB, (1) 
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where B and v are the large-scale (mean) magnetic and velocity fields cor- 
respondingly, a is the helicity coefficient and /? is the turbulent difFusivity. 

In the linear approximation the time dependent part can be factored out 
trivially: 

B(r,t) = B(r)e^*, (2) 

where 7 is an imaginary parameter to be found. 

We consider the aO-dynamo model and suppose that magnetic field gen- 
eration occurs in a thin spherical shell corresponding to the convective shell 
of the Sun. We also suppose that the differential rotation is more intensive 
than the mean helicity and that the radial gradient of the mean angular 
velocity G = weakly depends on latitude 6. 

To treat the problem, it is convenient to present the axisymmetric mag- 
netic field as a superposition of the toroidal and poloidal components and 
express the latter component as follows: Bp = Vx (0, 0, A) (we use the 
spherical system of coordinates and measure off the latitude from the equa- 
tor), where A a the vector potential. 

After averaging Eq.([l|) over the shell, one obtains 



7A = Q(e)B + ^ 



Z3z(^cos( 

cos a da 

1 d 



cos 9 dO 



{B cos 6) 



(3) 



(4) 



Note, that all quantities have been rescaled which leaded to the simple 
dimensionless form of the equations. Parameter D is the dimensionless dy- 
namo number which is supposed to be negative and numerically large. Func- 
tions a{6) and G{9) are measured in units of their maximum values and they 
are certainly unknown explicitly. Observations of the solar magnetic activ- 
ity and helioseismological data (see e.g. Schou et al., 1998) enable us to 
estimate them with a good accuracy {e.g., Belvedere et al., 2000). In what 
follows, we will suppose a{9) and G{9) to be reasonable smooth functions 
subject to the following symmetry constraints 

a{9) = -a{-9), and G{9) = G{-9). 

One of the advantages of the asymptotic method we apply is that the ex- 
plicit knowledge of the functions is not required and the final results can be 
expressed in a quite general form. To make explicit estimates, we use the 
following simple form of the functions: a{9) = sm9 and G = 1. 
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It is convenient to rewrite Eqs.(^, ^) in the following symbolic form: 



nh = -fh, (5) 
where we introduced a two-component complex function 



m = : : (6) 




and a matrix differential operator Ti. which is well-defined by Eqs.(^, ^) 

/ A ai9) 



n 



\ -DG{e)-^ cos A 



(7) 



where a notation for the azimuthal part of the Laplacian is introduced for 
brevity A = ^ ^ cos 6. Note that operator (|7|) is nonhermitian and, 

thus, its eigenvahies are complex in general. The adjoint operator has the 
following form 

/ A DGie)^ cose \ 
= . (8) 

V «(^) A J 

Note that operators Ti and are defined on the segment 6* G [— 7r/2, 7r/2]. 

The main idea of the WKB solution obtained in the above cited papers 
was to expand both the eigenvectors and eigenvalues into an asymptotic 
series in terms of parameter e = which was supposed to be small. 

Note, that magnetic field generation occurs only due to the helicity and dif- 
ferential rotation mechanisms [functions a{e) and G{6) cos{9)]. Remarkably, 
the two physically different mechanisms collapse into one universal function 
d{9) = a{9)G {e) cos{e) in the framework of the asymptotic solution. This 
function is relatively small near the equator and the poles. The magnetic 
field appears in these regions mainly due to the diffusion from middle lati- 
tudes where the generation is efficient. The solution was obtained under the 
assumption that the magnetic field generated in different hemispheres inde- 
pendently. This is true if the interplay between the southern and northern 
dynamo waves is negligible. 

In the framework of the WKB approach, we present the solution n in 
the northern hemisphere as a product of a fast oscillating exponent and a 
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iS{9) 



(9) 



slowly varying amplitude no: 

n{9) = no(6')exp 

The eigenvalues are seeking in the form of the asymptotic series 

-f = e-^To + e-^Ti,n, (10) 

where integer index n parameterizes the discrete eigenvalues. 

Using the WKB approximation technique, one can obtain the following 
explicit expression for the amplitude 



1 



cos ( 



ie~'^k{6) cos ( 



(11) 



where k{6) = dS{9)/d9 which satisfies the following Hamilton- Jacobi equa- 
tion: 



TQ + k'^{9) - ia{9)k{9) = 0, 



(12) 



where a = a{9)G{9) cos 9. This algebraic equation possesses four branches 
of roots. It is impossible to construct a smooth solution using any individual 
branch. Such a solution appears only by matching two different branches. A 
continuous crossover from one branch to the other is possible only for some 
unique values of the spectral parameter. This condition pins the value of 7 
and determines the desired spectrum: 



28/3 "max "= 



(13) 



and 



Ti^n = 3ik' {9') [n + 1/2], n = 0,1,2,..., (14) 

where 9' is the matching point which is the point of maximum of function 
a{9). The leading mode corresponds to the value n = 0. Function cr{9) in 
Eq.([lT|) has the following explicit form: 



CTnid) = exp 




l,n 



ik' [1 + 



2r 



To + k' 



2ik 



a 



d9 



(15) 



2(ro + A:^) 



Equations (|9| — |^) determine the asymptotic solution completely. We use 
these explicit formulae throughout the paper. 
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3 Equatorial region 



The crucial assumption of the WKB approximation is that the amplitude of 
the eigenvector varies in space much slower than the exponential factor [see 
Eq.Q]. This implies: 

dn^ 1 ^ 
dO e 

One can easily check that in the vicinity of the equator 0^1, the following 
estimate holds: 

Thus, in domain 9 < e, the applicability of the WKB approximation becomes 
questionable. 

Note, that at ^ = 0, two different branches k{6) merge. In the WKB 
theory, such a point is called "turning point" . It is more rule than exception 
that a WKB solution diverges at a turning point. Using explicit expressions 
— iiideed, observe that out solution diverges at ^ = 0. It is possible 



to show that 

a{e) « 0-1/4^ for 9^0. 
Moreover, one can see that k'{9) ~ 0-^/2 di^gr 

ges as well. 

Note, that these infinities have no physical meaning and appear due to 
unjustified approximations. The true solution is indeed a smooth function 
everywhere including the equator. This can be easily seen by expanding 
equations (^, ^) in the vicinity of 6* = 0: 

-fA{9) = a'{0)9B{9) + A"{9), (16) 

jB{9) = -DA' {9) + B"{9). (17) 

In these equations 7 plays the role of an independent external parameter. 
In the region 9 < e, the solution can be written as a superposition of two 
waves 

_ fA, + SA,i9)\^^e,(A2 + SA2{9)\^, 
-[B, + 6B,i9))'' ^[B2 + dB2{9))'' ' ^'"^^ 

where Ai^2 and Bi^2 are some constants to be determined by matching of 
( |T8| ) with the solution in the main domain. The ^-dependent corrections to 
the amplitudes in Eq. (lT8|) can be easily found explicitly. It is possible to 
show, that for 9^e these corrections are small and can be safely neglected 




in the leading approximation. Solution (18) is a perfectly smooth function 
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We denote the solution in the northern hemisphere as fi{9) and in the 
southern hemisphere as s{d). Let us note, that the first term in Eq.(|l8|) 
describes a plane wave propagating towards the equator, while the second 
one describes the wave propagating polewards and decaying exponentially. 
The first term can be matched with the WKB solution in the northern hemi- 
sphere n{6). The only way to match the second one is to suppose that the 
wave generated in the southern hemisphere s{0) penetrates to the northern 
hemisphere. In the vicinity of the equator the exponential factors in Eq.(p^) 
and in the WKB solution coincide. Equating the WKB eigenvector no with 
the first exponential term coefficient in ( [l8| ) and sq with the second one, we 
can express the coefficients Ai^2 and Bi^2 through the WKB parameters. 

Note that there are four different branches of roots k{9) of Eq.(^). Two 
of them describe the dynamo wave in the main domain. As it was shown by 
Galitski and Sokoloff (1999), the third branch describes the wave reflected 
from the pole. One can check that the fourth branch left is the only possible 
candidate to describe the wave propagating over the equator. Only this 
branch decays all the way down to ^ = —tt/2. 

Using the asymptotic expansions in the WKB solution for the case <C 1, 
we obtain: 

»=(,-.,e.»/..<V..j«P(-V7<') (19) 
and the adjoint solution in the southern hemisphere: 



^""^=1^ ae(-W3),(V4) )e^p[-VrO), (20) 

where 0i ~ e is the matching point and the following real constants have 
been introduced for brevity: 



a = 2^/33Vya'(0)a^/f. (21) 

and 

b = 3V22("V3) -1/3^ ^22) 

Now, we can estimate the decrement of spatial decay of the dynamo 
wave straying into the other hemisphere. Indeed, the characteristic latitude 
of its propagation ahead is l/^/Re'y. For the case D = —10^, or e = 0.1, we 
can estimate its leading order approximation using formulae ([To|) and ( |l3| ) 
as 0.28, or 16°. Notice that for D = —10^ this value is of order 7°. These 
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numbers look quite reasonable in view of the observational data obtained by 
Harvey (1992). She found the reversed polarity active regions that violate 
Hale's polarity law. Such regions appear in the final phase of the solar cycle 
mainly at low latitudes, near the solar equator. These results give signatures 
of the straying dynamo wave propagating from the counterpart hemisphere. 

4 Symmetry properties 

To formalize the subsequent calculations, let us introduce the following sym- 
metry operator P which we define as follows: 



right-hand side commutes with P. Thus, its eigenfunctions, i.e. magnetic 
field B, can be classified by the parity index p. Value p = +1 corresponds 
to the quadrupole solution, while value p = —1 corresponds to the dipole 
one. Eigenvalues 7 corresponding to different parities do not coincide. 

Since, Eqs.(|3[§) [or equivalently Eq.(P)] follow from Eq.(||), they inherent 
the symmetry properties of the initial equations. Let us note that these 
reduced equations involve a toroidal component of the physical magnetic 
field and a component of the gauge field. These two fields have opposite 
parities. Thus, when dealing with composite objects like (^) which contain 
the both fields, it is necessary to take into account this difference. Let us 
introduce the following unitary matrix: 



Pf{e) = fi-o). 



This operator has the following obvious eigenvalues (parity): 



p = ±l. 





(23) 



The dipole solution satisfies the following condition 



PUd = d 



(24) 



and consequently the quadrupole one is defined by 



PUq 



= -Q- 



(25) 
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The corresponding equations are written as 



= Jdd and Hq = q. (26) 

We also define the eigenvectors for the adjoint operator @ as 

n^cT = J* and n'^^ = 7*g". (27) 

Since the operator is nonhermitian, the eigenvectors for the mutually adjoint 
operators do not coincide. Let us note here that the eigenvalues may exist 
only in the form of complex conjugated pairs. Formally, this follows from 
the fact that operator IhL as well as the magnetic field are real. Below we are 
mainly interested in the structure of eigenfunctions. However, the splitting 
of the eigenvalues is of some interest as well. 

To find the splitting of the eigenvalues, we make use of a method used to 
solve a well-known quantum-mechanical problem. Namely, we observe that 
the problem at hand is very similar to the double-well problem in quantum 
mechanics. In the latter, a particle in a symmetric one-dimensional potential 
is studied. The potential consists of two quantum wells separated by a high 
barrier. If one neglects the possibility of the penetration through the barrier 
the eigenvalues are degenerated and they can be calculated in a well (say in 
the right one) with the help of the WKB approximation. The eigenfunctions 
(wave functions) decay exponentially far from the well. If one takes into ac- 
count a finite probability of the barrier penetration, the degeneracy is lifted 
and the energy levels split into pairs corresponding to the symmetric and 
antisymmetric solutions. The quantity of interest is the energy difference 
between the two lowest eigenstates, which corresponds to the tunneling rate 
through the double well barrier. The quantity is often small and difficult to 
calculate numerically, especially when the potential barrier between the two 
wells is large. However, using the WKB eigenfunctions obtained for each 
quantum wells, it is easy to construct approximate symmetric and antisym- 
metric solutions explicitly. The subsequent calculation of the tunneling rate 
is straightforward and simple (see e.g. Landau and Lifshitz, 1968). 

In our problem we are dealing with the southern and northern domains of 
generation separated by the equatorial region in which magnetic field gener- 
ation is weak. This equatorial region corresponds to the barrier. Neglecting 
the interaction of the dynamo waves generated in different hemispheres, one 
can obtain the WKB eigenfunctions and eigenvalues explicitly [see Eqs.(^ — 



15)]. Using the solution in the northern hemisphere n{9) and in the southern 



10 



one s{9), we follow the classical Lifshitz solution of the quantum problem 
and construct the dipole and quadrupole eigenfunctions as follows 

d{0) = ±[ni9) + m] (28) 

and 

m = ^ im - m] , (29) 

where factor l/\/2 is introduced in order to preserve the norms of the eigen- 
vectors. 

Let us note, that the solution in the southern hemisphere can be obtained 
easily by writing 

s{e) = PUn{e). 

To proceed further, let us introduce the following "inner product" of two 
vector functions / and g: 

7r/2 

(/,5)= / [h{0)9m + f2{O)gm]^osede, (30) 



where 



By multiplying equation Tin = 7on on and Eq.(^) on n, we obtain 

7d-lo = ^ ^ .J ; ^ (31) 

and 

Let us note here that the dynamo- wave generated in the southern hemisphere 
s{9), if exists, is exponentially small in the northern one. Thus, we conclude: 

|(s-«*,n)|«|(n-,n)|. 

Hence 



(tf * , nj ~ (q^* , n) ~ {n"* , n) 
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With a good accuracy, we can neglect the corresponding difference in the de- 



nominators of expressions ( |3l] ) and (32) and obtain the following important 
identity: 

((n^^y , n) - (^,nn) 

This equation brings the dynamo problem into direct correspondence with 
the quantum one. 

Using the explicit expressions for the operators "ti and 7^^, we obtain the 
following relation: 



^(n-,n)A7* 



7r/2 



dQ cos Q 



As\nx + (cos Q sf) ni + As%n2 
du 



-si (Ani) + L»4— (cos^m) -s?(An2) 
du 



(34) 



where indexes "1" and "2" correspond to the upper and lower components 
of a vector and s"" means the solution of the adjoint equation in the south- 



ern hemisphere. Integral Eq.(34) can be easily evaluated by parts and the 
splitting is expressed as boundary terms. 



1 



dsl 



ds9, 



,dni 



. dn-) 



-in"* , n) A7* = ni— f + n2-^ - s?— ^ - s^-f + Dms^ 



d9 



dO 



de 



dO 



(35) 



5 Eigenvalue splitting 

In this section we calculate the eigenvalue splitting explicitly using the WKB 
solution obtained previously (see Sec. 2). First of all, we find the eigenso- 
lution of the adjoint operator (^). 

Let vector n{9), with the upper component ni{9) and lower n2{6), be 
the WKB solution of Eq.(||) [see Eqs. (^ — 15)]. Then, the corresponding 
solution of the adjoint equation can be obtained by the following replace- 
ments: 

"-1,2 "-2,1, 

Fo^rS and Fi^F^. 

Function k{9) satisfying the Hamilton-Jakobi equation ([l^ ) should be re- 
placed by the following value: 



kie) -k*{e). 
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After the appropriate replacements are made, the adjoint solution takes the 
form: 

' -ie~^k* (9) cos e\ a* (6) [ iS* (6) " 

k*{ey ) cos (9 



exp 



(36) 



At this point, we can calculate the eigenvalue splitting using the explicit 
expressions obtained above. We start with the evaluation of the inner prod- 
uct (n" , n). Let us note, that in the double- well problem the corresponding 
product is nothing but the norm of the wave-function and, thus, is equal to 



unity. Using Eqs.(ll) and ( pq ) we get: 

tt/2 



{n" , n) 



1 



'- (fo + k^{e)) ik*{0) + ik{e) (to + k^e))* 

exp J y" Im k{e')de' I . 



(37) 



Since the integrand in the above expression contains a Gaussian exponent 
with a sharp maximum, it is possible to evaluate the integral by the Laplace 
method. In the saddle point approximation, we obtain 



4Imr 







vre 



Im k' 



k{e^) exp I J lmk{e)d9 



(38) 



where 9^, is the point where action S{9) reaches its maximum, i.e. the point 
at which ImA;(^=K) = 0. This point has already been found explicitly by 
Kuzanyan and Sokoloff (1995) as follows 



9\/3 



0.8052. 



16^/2V^/3- 1 



(39) 



Let us note that quantity ( |38D is a real, positive and exponentially large 
number. Thus, the eigenvalue splitting is exponentially small which is in 
accord with the familiar result for the double well problem. 

Using expressions (p^)~(p2|) , we obtain the following results for boundary 
terms (ISSl): 



ni 



dsf 



,dn2 
' d9 



2ab, 



-Re 



(40) 
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n2- 



de 



2ab^ 



-Re 



.(-W6), 



(41) 



DniS2 



'1 (- 



■j7r/3) 



(42) 



=0 

The two first boundary terms give real positive contributions to the spht- 
ting. This is not very surprising, since the corresponding terms come from 
the hermitian part of operator Ti.. What is more remarkable is that the 
matching point 6i drops out of the final result for these terms. The third 
"nonhermitian" term gives a nonvanishing contribution to the imaginary 
part of the eigenvalue difference and explicitly depends upon the matching 
point. Since 9i ~ e, we see that this contribution is parametrically small 
compared to the first two ones. 

we obtain the following result: 



Using Eqs.(gg, |l 



(n" , n) Re A7 



1 39/^ 
£3 21/3 



a 



4/3 



d'(0) 



and 



(n" , n) Im A7 



219/63I/ 




(43) 



(44) 



where norm (n"* , n) was calculated in Sec. 4 [see Eq. ([3^)]. The matching 
point 61 can be calculated self-consistently as a point at which the phase- 
shifts for different asymptotics coincide [cf. Galitski and Sokoloff (1999)]. 

As we have seen above, the WKB approximation breaks down near the 
equator. This happens because 9 = is the turning point for our solution. 
Thus, straightforward evaluation of (^5|) using the WKB formulae is not 
possible. Note that in the usual quantum problem such a difficulty does not 
arise since the boundary is located far from the turning point. The WKB 
solution can be applied directly and the exponential term coefficient can be 
found easily. However, it is important to realize that the method itself is 
limited by the exponential accuracy. The numerical part of the exponen- 
tial term coefficient should not be trusted even if found. There are other 
asymptotical methods allowing the calculation of the tunneling rate, which 
give slightly different results. Indeed, the exponential factor is the same for 
all these methods. For a detailed discussion of the asymptotical methods in 
quantum double-well problems see e.g. Cooper et al. (1995). However, as 
we have already mentioned, the method used does not allow us to obtain 
the correct value of the exponential term coefficient. We conclude that the 
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accuracy of Eqs.(43, 44) already exceeds the accuracy of the method. Hence, 
the exact value of ~ e is not important. 

Bearing these observations in mind, let us write the final result in the 
following form: 



Re A7 

I70I 



\/e exp < - 

e Jo 



\lm k{9)\de 



(45) 



and 



Im A7 

I70I 



-eexp 



2 

e Jo 



\lmk{e)\d9\ 



(46) 



6 Discussion 



The overall picture of the magnetic field in the Parker model can be now 
summarized as follows. A strong dynamo wave appears at middle latitudes. 
The point where the generation sources reach maximum is shifted from the 
point where the helicity coefficient is maximal. This happens because the 
differential rotation and helicity mechanisms are both important on equal 
footing. In the framework of the asymptotic analysis, the two mechanisms 
manifest themselves through a universal function. The asymptotic WKB 
theory allows to derive explicit expressions describing the generated dynamo 
waves. It appears that at middle latitudes the main wave propagates equa- 
torwards. At very high latitudes the wave changes its direction of motion 
and propagates towards the pole and reflects from it. In the vicinity of the 
equator there is a transition zone ~ e where the WKB solution becomes 
progressively worse and the crossover from the asymptotic behaviour (^ to 
( p!8| ) occurs. The wave does not vanish completely at the equator but rolls 
over it propagating in the southern part. Beyond the southern transition 
zone 6 ~ — e, the WKB solution becomes applicable again and it describes 
a very weak decaying dynamo wave. 

Let us mention here the observational results which refer to the time 
of the Maunder minimum (see Sokoloff and Nesme-Ribes, 1994; Ribes and 
Nesme-Ribes, 1993). During the Maunder minimum, the solar magnetic 
activity was very weak. At the end of the minimum, some signs of the ac- 
tivity appeared in the southern hemisphere. The southern wave of magnetic 
activity looked normal and propagated equatorwards. In the northern part 
there was almost no activity at all, except an unusual wave which existed 
only near the equator and propagated away from it. 
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Certainly, our simple linear asymptotic theory can not suggest any expla- 
nations for the Maunder minimum itself. However, we think that the weak 
equatorial wave detected at the epoch of the minimum corresponds to the 
equatorial wave which appears naturally in the framework of our asymptotic 
analysis. On the modern Sun, this equatorial wave, if exists, is screened by 
the background of the main dynamo wave. There axe observational signa- 
tures of some magnetic field tracers such as extreme ultraviolet (EUV) lines 
to propagate polewards near the solar equator {e.g., Benevolenskaya et al., 
2001). They may indicate the presence of such a straying wave propagating 
from the other hemisphere. Notice, that the analysis of the observational 
data by Harvey (1992) revealed that the reversed polarity active regions that 
violate the Hale's polarity law are located mainly at low latitudes. Further 
observations of the equatorial region would be highly desirable to clarify this 
question. 

We have shown that a small interaction between the northern and south- 
ern waves yields a splitting of the eigenvalues corresponding to the dipole 
and quadrupole configurations. It appears that the growth rate correspond- 
ing to the dipole configuration exceeds the quadrupole one. This fact gives 
a very tentative indication that the dipole configuration is to dominate the 
quadrupole one after a sufficiently long time. Let us estimate this time using 
the following trial parameters: a{0) = sin 9 and D = —10^. In this case 

^ ^ 0.03 

l7o| 

If we suppose that (27r/Im7o) ~ 22yr, where Im7o = |7o|'\/3/2. Then we 
have the following estimate for the time of dominance of the dipole mode 
over the quadrupole one r = 1/Re A7 ~ 100 yr. This is comparable with the 
time of recovery of the solar magnetic field generation (Ribes and Nesme- 
Ribes, 1993) from the Maunder Minumum when some activity was confined 
to one southern hemisphere (equatorial asymmetry). 

The difference in the imaginary part of the eigenvalues may give rise to 
an additional weak period of oscillations of the magnetic field in the non- 
linear regime. The corresponding estimate for the period of such a global 
modulation is T2 = 27r/Im A7 ~ 2 • 10'^ yr. 

Let us stress that the applicability of the simple Parkler's migratory dy- 
namo model is limited. A two-dimensional and non- linear generalizations 
are required to provide a more accurate description. However, some simple 
estimates can be done now. Namely, using a two-dimensional asymptotic 
solution, e.g., obtained recently by Belvedere et al. (2000), we can estimate 
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the overlap of the dynamo waves generated in different hemispheres. Be- 
cause the r.h.s. of equation (34) then turns to a 2D integral, we expect 
that the corresponding overlap exceeds significantly that one in the one- 
dimensional model. This may result in a larger splitting of the eigenvalues. 
The corresponding period of secondary oscillations T2 should increase. This 
qualitative result suggests a possible explanation for the well-known Gleiss- 
berg cycle. 
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